The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.

Biochemical changes in the frontal lobe of HIV-infected individuals detected by magnetic resonance spectroscopy

We have developed a proton magnetic resonance spectroscopy method that selectively can sample cortical gray matter and adjacent white matter in the frontal lobe. We have used this approach to study a group of patients (n = 7) infected with HIV and clinical manifestations of the AIDS dementia complex (ADC), a group of patients (n = 8) infected with HIV without any indications of ADC, and seven controls. The patients without ADC had a statistically significant increase in the ratio of myo-inositol to creatine in white matter compared with normal controls. In contrast, the group of patients with ADC had almost normal levels of myo-inositol to creatine in both gray matter and white matter and showed a statistically significant decrease in the N-acetylaspartate to creatine ratio in gray matter compared with either the normal controls or the patients without ADC. Patterns of spectral abnormalities correlated with neuropsychological measures of frontal lobe dysfunction, suggesting that the evaluation of frontal lobe metabolism by magnetic resonance spectroscopy can play a role in the early detection of ADC, in determining its progression, and in assessing responses to therapeutic interventions.

Global variation in the genetic and biochemical basis of diamondback moth resistance to Bacillus thuringiensis

Insecticidal proteins from the soil bacterium Bacillus thuringiensis (Bt) are becoming a cornerstone of ecologically sound pest management. However, if pests quickly adapt, the benefits of environmentally benign Bt toxins in sprays and genetically engineered crops will be short-lived. The diamondback moth (Plutella xylostella) is the first insect to evolve resistance to Bt in open-field populations. Here we report that populations from Hawaii and Pennsylvania share a genetic locus at which a recessive mutation associated with reduced toxin binding confers extremely high resistance to four Bt toxins. In contrast, resistance in a population from the Philippines shows multilocus control, a narrower spectrum, and for some Bt toxins, inheritance that is not recessive and not associated with reduced binding. The observed variation in the genetic and biochemical basis of resistance to Bt, which is unlike patterns documented for some synthetic insecticides, profoundly affects the choice of strategies for combating resistance.

Toward an understanding of the biochemical and pharmacological complexity of the saliva of a hematophagous sand fly Lutzomyia longipalpis

The saliva of blood-sucking arthropods contains powerful pharmacologically active substances and may be a vaccine target against some vector-borne diseases. Subtractive cloning combined with biochemical approaches was used to discover activities in the salivary glands of the hematophagous fly Lutzomyia longipalpis. Sequences of nine full-length cDNA clones were obtained, five of which are possibly associated with blood-meal acquisition, each having cDNA similarity to: (i) the bed bug Cimex lectularius apyrase, (ii) a 5′-nucleotidase/phosphodiesterase, (iii) a hyaluronidase, (iv) a protein containing a carbohydrate-recognition domain (CRD), and (v) a RGD-containing peptide with no significant matches to known proteins in the blast databases. Following these findings, we observed that the salivary apyrase activity of L. longipalpis is indeed similar to that of Cimex apyrase in its metal requirements. The predicted isoelectric point of the putative apyrase matches the value found for Lutzomyia salivary apyrase. A 5′-nucleotidase, as well as hyaluronidase activity, was found in the salivary glands, and the CRD-containing cDNA matches the N-terminal sequence of the HPLC-purified salivary anticlotting protein. A cDNA similar to α-amylase was discovered and salivary enzymatic activity demonstrated for the first time in a blood-sucking arthropod. Full-length clones were also found coding for three proteins of unknown function matching, respectively, the N-terminal sequence of an abundant salivary protein, having similarity to the CAP superfamily of proteins and the Drosophila yellow protein. Finally, two partial sequences are reported that match possible housekeeping genes. Subtractive cloning will considerably enhance efforts to unravel the salivary pharmacopeia of blood-sucking arthropods.

Interaction of Fibroblast Growth Factor-2 (FGF-2) with Free Gangliosides: Biochemical Characterization and Biological Consequences in Endothelial Cell Cultures

Exogenous gangliosides affect the angiogenic activity of fibroblast growth factor-2 (FGF-2), but their mechanism of action has not been elucidated. Here, a possible direct interaction of sialo-glycolipids with FGF-2 has been investigated. Size exclusion chromatography demonstrates that native, but not heat-denatured, 125I-FGF-2 binds to micelles formed by gangliosides GT1b, GD1b, or GM1. Also, gangliosides protect native FGF-2 from trypsin digestion at micromolar concentrations, the order of relative potency being GT1b > GD1b > GM1 = GM2 = sulfatide > GM3 = galactosyl-ceramide, whereas asialo-GM1, neuraminic acid, and N-acetylneuramin-lactose were ineffective. Scatchard plot analysis of the binding data of fluorochrome-labeled GM1 to immobilized FGF-2 indicates that FGF–2/GM1 interaction occurs with a Kd equal to 6 μM. This interaction is inhibited by the sialic acid-binding peptide mastoparan and by the synthetic fragments FGF-2(112–129) and, to a lesser extent, FGF-2(130–155), whereas peptides FGF-2(10–33), FGF-2(39–59), FGF-2(86–96), and the basic peptide HIV-1 Tat(41–60) were ineffective. These data identify the COOH terminus of FGF-2 as a putative ganglioside-binding region. Exogenous gangliosides inhibit the binding of 125I-FGF-2 to high-affinity tyrosine-kinase FGF-receptors (FGFRs) of endothelial GM 7373 cells at micromolar concentrations. The order of relative potency was GT1b > GD1b > GM1 > sulfatide a = sialo-GM1. Accordingly, GT1b,GD1b, GM1, and GM2, but not GM3 and asialo-GM1, prevent the binding of 125I-FGF-2 to a soluble, recombinant form of extracellular FGFR-1. Conversely, the soluble receptor and free heparin inhibit the interaction of fluorochrome-labeled GM1 to immobilized FGF-2. In agreement with their FGFR antagonist activity, free gangliosides inhibit the mitogenic activity exerted by FGF-2 on endothelial cells in the same range of concentrations. Also in this case, GT1b was the most effective among the gangliosides tested while asialo-GM1, neuraminic acid, N-acetylneuramin-lactose, galactosyl-ceramide, and sulfatide were ineffective. In conclusion, the data demonstrate the capacity of exogenous gangliosides to interact with FGF-2. This interaction involves the COOH terminus of the FGF-2 molecule and depends on the structure of the oligosaccharide chain and on the presence of sialic acid residue(s) in the ganglioside molecule. Exogenous gangliosides act as FGF-2 antagonists when added to endothelial cell cultures. Since gangliosides are extensively shed by tumor cells and reach elevated levels in the serum of tumor-bearing patients, our data suggest that exogenous gangliosides may affect endothelial cell function by a direct interaction with FGF-2, thus modulating tumor neovascularization.

Functional and Biochemical Analysis of the C2 Domains of Synaptotagmin IV

Synaptotagmins (Syts) are a family of vesicle proteins that have been implicated in both regulated neurosecretion and general membrane trafficking. Calcium-dependent interactions mediated through their C2 domains are proposed to contribute to the mechanism by which Syts trigger calcium-dependent neurotransmitter release. Syt IV is a novel member of the Syt family that is induced by cell depolarization and has a rapid rate of synthesis and a short half-life. Moreover, the C2A domain of Syt IV does not bind calcium. We have examined the biochemical and functional properties of the C2 domains of Syt IV. Consistent with its non–calcium binding properties, the C2A domain of Syt IV binds syntaxin isoforms in a calcium-independent manner. In neuroendocrine pheochromocytoma (PC12) cells, Syt IV colocalizes with Syt I in the tips of the neurites. Microinjection of the C2A domain reveals that calcium-independent interactions mediated through this domain of Syt IV inhibit calcium-mediated neurotransmitter release from PC12 cells. Conversely, the C2B domain of Syt IV contains calcium binding properties, which permit homo-oligomerization as well as hetero-oligomerization with Syt I. Our observation that different combinatorial interactions exist between Syt and syntaxin isoforms, coupled with the calcium stimulated hetero-oligomerization of Syt isoforms, suggests that the secretory machinery contains a vast repertoire of biochemical properties for sensing calcium and regulating neurotransmitter release accordingly.

Phosphorylation Controls Timing of Cdc6p Destruction: A Biochemical Analysis

The replication initiation protein Cdc6p forms a tight complex with Cdc28p, specifically with forms of the kinase that are competent to promote replication initiation. We now show that potential sites of Cdc28 phosphorylation in Cdc6p are required for the regulated destruction of Cdc6p that has been shown to occur during the Saccharomyces cerevisiae cell cycle. Analysis of Cdc6p phosphorylation site mutants and of the requirement for Cdc28p in an in vitro ubiquitination system suggests that targeting of Cdc6p for degradation is more complex than previously proposed. First, phosphorylation of N-terminal sites targets Cdc6p for polyubiquitination probably, as expected, through promoting interaction with Cdc4p, an F box protein involved in substrate recognition by the Skp1-Cdc53-F-box protein (SCF) ubiquitin ligase. However, in addition, mutation of a single, C-terminal site stabilizes Cdc6p in G2 phase cells without affecting substrate recognition by SCF in vitro, demonstrating a second and novel requirement for specific phosphorylation in degradation of Cdc6p. SCF-Cdc4p– and N-terminal phosphorylation site–dependent ubiquitination appears to be mediated preferentially by Clbp/Cdc28p complexes rather than by Clnp/Cdc28ps, suggesting a way in which phosphorylation of Cdc6p might control the timing of its degradation at then end of G1 phase of the cell cycle. The stable cdc6 mutants show no apparent replication defects in wild-type strains. However, stabilization through mutation of three N-terminal phosphorylation sites or of the single C-terminal phosphorylation site leads to dominant lethality when combined with certain mutations in the anaphase-promoting complex.

Genetic and Biochemical Characterization of the Yeast Spo12 Protein

We have performed a genetic and biochemical analysis of the SPO12 gene, which regulates meiotic nuclear divisions in budding yeast. When sporulated, spo12 mutants undergo a single meiotic nuclear division most closely resembling meiosis II. We observed that Spo12 protein is localized to the nucleus during both meiotic divisions and that Clb1-Cdc28, Clb3-Cdc28, Clb4-Cdc28, and Clb5-Cdc28 kinase activities during meiosis were not affected by a spo12 mutation. Using two-hybrid analysis, we identified several genes, three of which are meiotically induced, that may code for proteins that interact with Spo12p. We also observed that two genes, BNS1 (Bypasses Need for Spo12p), which has homology to SPO12, and SPO13, whose mutant phenotype is like that of spo12, can partially suppress the meiotic defect of spo12 mutants when overexpressed. We found that Spo12p is also localized to the nucleus in vegetative cells and that its level peaks during G2/M. We observed that a spo12 mutation is synthetically lethal in vegetative cells with a mutation in HCT1, a gene necessary for cells to exit mitosis, suggesting that Spo12p may have a role in exit from mitosis.

Ocular cicatricial pemphigoid antigen: Partial sequence and biochemical characterization

Ocular cicatricial pemphigoid (OCP) is an autoimmune disease that affects mainly conjunctiva and other squamous epithelia. OCP is histologically characterized by a separation of the epithelium from underlying tissues within the basement membrane zone. Immunopathological studies demonstrate the deposition of anti-basement membrane zone autoantibodies in vivo. Purified IgG from sera of patients with active OCP identified a cDNA clone from a human keratinocyte cDNA library that had complete homology with the cytoplasmic domain of β4-integrin. The sera recognized a 205-kDa protein in human epidermal, human conjunctiva, and tumor cell lysates that was identified as β4-integrin by its reaction with polyclonal and monoclonal antibodies to human β4-integrin. Sera from patients with bullous pemphigoid, pemphigus vulgaris, and cicatricial pemphigoid-like diseases did not recognize the 205-kDa protein, indicating the specificity of the binding. These data strongly implicate a role for human β4-integrin in the pathogenesis of OCP. It should be emphasized that multiple antigens in the basement membrane zone of squamous epithelia may serve as targets for a wide spectrum of autoantibodies observed in vesiculobullous diseases. Molecular definition of these autoantigens will facilitate the classification and characterization of subsets of cicatricial pemphigoid and help distinguishing them from bullous pemphigoid. This study highlights the function and importance of β4-integrin in maintaining the attachment of epithelial cells to the basement membrane.

Manipulating the genetic identity and biochemical surface properties of individual cells with electric-field-induced fusion

A method for cell–cell and cell–liposome fusion at the single-cell level is described. Individual cells or liposomes were first selected and manipulated either by optical trapping or by adhesion to a micromanipulator-controlled ultramicroelectrode. Spatially selective fusion of the cell–cell or cell–liposome pair was achieved by the application of a highly focused electric field through a pair of 5-μm o.d. carbon-fiber ultramicroelectrodes. The ability to fuse together single cells opens new possibilities in the manipulation of the genetic and cellular makeup of individual cells in a controlled manner. In the study of cellular networks, for example, the alteration of the biochemical identity of a selected cell can have a profound effect on the behavior of the entire network. Fusion of a single liposome with a target cell allows the introduction of the liposomal content into the cell interior as well as the addition of lipids and membrane proteins onto the cell surface. This cell–liposome fusion represents an approach to the manipulation of the cytoplasmic contents and surface properties of single cells. As an example, we have introduced a membrane protein (γ-glutamyltransferase) reconstituted in liposomes into the cell plasma membrane.

Altering the biochemical state of individual cultured cells and organelles with ultramicroelectrodes

We describe an efficient technique for the selective chemical and biological manipulation of the contents of individual cells. This technique is based on the electric-field-induced permeabilization (electroporation) in biological membranes using a low-voltage pulse generator and microelectrodes. A spatially highly focused electric field allows introduction of polar cell-impermeant solutes such as fluorescent dyes, fluorogenic reagents, and DNA into single cells. The high spatial resolution of the technique allows for design of, for example, cellular network constructions in which cells in close contact with each other can be made to possess different biochemical, biophysical, and morphological properties. Fluorescein, and fluo-3 (a calcium-sensitive fluorophore), are electroporated into the soma of cultured single progenitor cells derived from adult rat hippocampus. Fluo-3 also is introduced into individual submicrometer diameter processes of thapsigargin-treated progenitor cells, and a plasmid vector cDNA construct (pRAY 1), expressing the green fluorescent protein, is electroporated into cultured single COS 7 cells. At high electric field strengths, observations of dye-transfer into organelles are proposed.

Biochemical evolution II: Origin of life in tubular microstructures on weathered feldspar surfaces

Mineral surfaces were important during the emergence of life on Earth because the assembly of the necessary complex biomolecules by random collisions in dilute aqueous solutions is implausible. Most silicate mineral surfaces are hydrophilic and organophobic and unsuitable for catalytic reactions, but some silica-rich surfaces of partly dealuminated feldspars and zeolites are organophilic and potentially catalytic. Weathered alkali feldspar crystals from granitic rocks at Shap, north west England, contain abundant tubular etch pits, typically 0.4–0.6 μm wide, forming an orthogonal honeycomb network in a surface zone 50 μm thick, with 2–3 × 106 intersections per mm2 of crystal surface. Surviving metamorphic rocks demonstrate that granites and acidic surface water were present on the Earth’s surface by ∼3.8 Ga. By analogy with Shap granite, honeycombed feldspar has considerable potential as a natural catalytic surface for the start of biochemical evolution. Biomolecules should have become available by catalysis of amino acids, etc. The honeycomb would have provided access to various mineral inclusions in the feldspar, particularly apatite and oxides, which contain phosphorus and transition metals necessary for energetic life. The organized environment would have protected complex molecules from dispersion into dilute solutions, from hydrolysis, and from UV radiation. Sub-micrometer tubes in the honeycomb might have acted as rudimentary cell walls for proto-organisms, which ultimately evolved a lipid lid giving further shelter from the hostile outside environment. A lid would finally have become a complete cell wall permitting detachment and flotation in primordial “soup.” Etch features on weathered alkali feldspar from Shap match the shape of overlying soil bacteria.

Structure-based assignment of the biochemical function of a hypothetical protein: A test case of structural genomics

Many small bacterial, archaebacterial, and eukaryotic genomes have been sequenced, and the larger eukaryotic genomes are predicted to be completely sequenced within the next decade. In all genomes sequenced to date, a large portion of these organisms’ predicted protein coding regions encode polypeptides of unknown biochemical, biophysical, and/or cellular functions. Three-dimensional structures of these proteins may suggest biochemical or biophysical functions. Here we report the crystal structure of one such protein, MJ0577, from a hyperthermophile, Methanococcus jannaschii, at 1.7-Å resolution. The structure contains a bound ATP, suggesting MJ0577 is an ATPase or an ATP-mediated molecular switch, which we confirm by biochemical experiments. Furthermore, the structure reveals different ATP binding motifs that are shared among many homologous hypothetical proteins in this family. This result indicates that structure-based assignment of molecular function is a viable approach for the large-scale biochemical assignment of proteins and for discovering new motifs, a basic premise of structural genomics.

The biochemical properties of the ATPase activity of a 70-kDa heat shock protein (Hsp70) are governed by the C-terminal domains

The cytosolic 70-kDa heat shock proteins (Hsp70s), Ssa and Ssb, of Saccharomyces cerevisiae are functionally distinct. Here we report that the ATPase activities of these two classes of Hsp70s exhibit different kinetic properties. The Ssa ATPase has properties similar to those of other Hsp70s studied, such as DnaK and Hsc70. Ssb, however, has an unusually low steady-state affinity for ATP but a higher maximal velocity. In addition, the ATPase activity of Hsp70s, like that of Ssa1, depends on the addition of K+ whereas Ssb activity does not. Suprisingly, the isolated 44-kDa ATPase domain of Ssb has a Km and Vmax for ATP hydrolysis similar to those of Ssa, rather than those of full length Ssb. Analysis of Ssa/Ssb fusion proteins demonstrates that the Ssb peptide-binding domain fused to the Ssa ATPase domain generates an ATPase of relatively high activity and low steady-state affinity for ATP similar to that of native Ssb. Therefore, at least some of the biochemical differences between the ATPases of these two classes of Hsp70s are not intrinsic to the ATPase domain itself. The differential influence of the peptide-binding domain on the ATPase domain may, in part, explain the functional uniqueness of these two classes of Hsp70s.